Environmental Research Letters

AimHere we aim to disentangle species-specific bioclimatic drivers of forest site productivity and project their future dynamics, providing a spatially explicit basis for anticipating climate-driven shifts in productivity and their implications for forest carbon sequestration. LocationIberian Peninsula. Time period1985-2014 (calibration); 2071-2100 (projected under CMIP6 scenarios). Major taxa studied21 Iberian tree species. MethodsWe used Site Form (SF) maps derived from the Third Spanish National Forest Inventory, spatially interpolating plot-level SF estimates as a continuous productivity index and relating them to 25 bioclimatic variables. Multiple linear regression models were selected via complementary stepwise and subset regression and validated on independent hold-out data (80%/20% split). ResultsValidated [Formula] ranged from 0.46 (Quercus faginea) to 0.97 (Pinus pinaster); 17 of 21 species reached [Formula]. BI013 precipitation of the wettest month), not BI014, was the most frequently retained predictor (15/17); BI014 was retained in only (11/17 models with a near-even sign split. Combining projected changes in mean productivity and habitat extent under SSP5-8.5, fifteen of sixteen applicable species lose total productivity by 2071-2100, six -- including Fagus sylvatica and Betula alba -- collapsing to below 1% of their reference-period value; only Pinus pinaster gains, and only under the lowest-emission pathway (up to 175%) -- under SSP5-8.5 it too loses productivity, albeit less than any other species (35% of its reference-period value retained). Limiting warming to SSP1-2.6 spares Mediterranean pine and oak species but not Euro-Siberian and montane ones. Main conclusionsThese validated, extrapolation-aware models reveal a near-universal, climate-driven collapse in Iberian forest site productivity, with direct implications for the carbon-sink potential currently attributed to these forest types, and provide a route to dynamic, climate-aware carbon-uptake estimates for the region.

Businesses are increasingly encouraged to disclose their nature-related dependencies, impacts, risks and opportunities. A common component of sustainability reporting is screening operational sites for ecologically sensitivity to identify locations for further evaluation and action. However, with 600+ biodiversity metrics available, selecting and interpreting appropriate metrics remains challenging for business. We developed a simple screening framework informed by the Taskforce for Nature-Related Financial Disclosures guidance, grouping eleven widely used global biodiversity metrics into four complementary [&prime]baskets[&prime], representing different aspects of biodiversity. We created hypothetical but realistic mining, onshore wind energy and agricultural companies, to assess how metric choice, buffer size, scoring approach and sensitivity thresholds influence screening outcomes. Our basket framework consistently identified similar high-priority sites across metric combinations, but site rankings varied with methodological choices. We recommend clearer guidance on metric selection and application, alongside greater transparency from business regarding assumptions, methods and limitations when screening sites for ecological sensitivity.

Coastal forests are increasingly threatened by saturated soil and elevated salinity levels resulting from sea level rise, saltwater intrusion, and storm surges. In response to rising salinization and flooding, healthy coastal forests that rely on freshwater (both wetland forests and low-elevation upland forests) are transitioning into landscapes dominated by dead or dying trees, known as ghost forests. Situated among salt-tolerant shrubs and grasses, ghost forests eventually become marshes or open water. Here, our main objective was to quantify the dynamics and pathways of these forest landscape conversions, as well as the factors contributing to the changes, which is vital for understanding the progression of coastal ecosystem degradation and forecasting future changes. We focused first on identifying the best method to track forest landscape change by exploring the role of multiple remote sensing indices (i.e., multispectral, bi-seasonal, topographical, and phenological metrics) in enhancing the performance of deep learning models (convolutional neural networks, CNNs) for land cover classification in the coastal plain of North Carolina using surface reflectance of Landsat 8 and Sentinel-2 images. Then, we used the best available data (Landsat 8) to understand long-term change and identify patterns of land cover change from 1985 to 2021. Our study reveals that incorporating phenology and topographical indices enhances the separability of the ghost forests class from all other vegetation classes. In our assessment, the higher-resolution Sentinel-2 data (F1 Score = 96.3) outperformed Landsat images (F1 score = 93.4) for the 2021 co-available year. However, Landsat remains an important tool used due to its long-term data record. Therefore, we used Landsat to determine that 21% of forests were lost between 1985 and 2021, and that the rate of loss is increasing. Between 2010 and 2021, 23,876 ha of forest were converted to marsh, ghost forest, and shrub, which is 1.5 times higher than the 16,968 ha lost between 1985 and 2010. These conversions from forest to ghost forest and marshes were driven primarily by proximity to the channel, salinity, and the increasing rate of relative sea level rise (RSLR), which are the key environmental drivers of observed changes. By quantifying these changes, we highlight regions most vulnerable to environmental stressors, providing a basis for targeted conservation strategies.

Offshore wind farms are rapidly expanding to meet growing demands for renewable energy, with development expected to extend further offshore into deeper waters. This expansion requires a robust understanding of the long-term ecological consequences of offshore wind farms (OWFs) and how these may interact with ongoing climate change. We used the coupled hydrodynamic-ecosystem-biogeochemical water-column model (GOTM-ERSEM-BFM) to investigate ecosystem-wide responses to environmental changes associated with OWFs and climate warming. Specifically, we examined OWF-related scenarios of reduced benthic suspension-feeding activity, representing potential effects of contaminant emissions from OWFs, and reduced wind forcing, together with increased sea surface temperature. The scenarios were simulated individually and in combination to explore potential interactive effects. These scenarios were simulated at two contrasting locations in the North Sea, representing a well-mixed coastal site and a seasonally stratified offshore site. The coastal site exhibited comparatively modest ecosystem responses across the scenarios, whereas responses were generally stronger at the deeper offshore site. At the offshore site, changes in stratification altered vertical nutrient dynamics and contributed to pronounced differences in ecosystem responses between the surface and bottom layers. Our results demonstrate that ecosystem responses to OWF-related and climate-driven environmental changes are strongly dependent on local environmental conditions, suggesting that ecological consequences may differ substantially as wind farm development expands into deeper offshore environments.

The global market concentrates extractive pressure on lands held by Indigenous peoples, including peoples living in isolation, for whom free, prior and informed consent cannot be obtained and protection must therefore rest on territorial instruments. Halmahera, Indonesia, holds some of the worlds largest lateritic nickel reserves beneath a lowland rainforest inhabited by the Hongana Manyawa, yet the trajectory of land use and cover change across the island has not been quantified. We characterised land use and cover change over the 17,437 km2 island between 2014 and 2024 using MapBiomas time series, and projected a business-as-usual scenario to 2054 with a stochastic cellular-automata model implemented in Dinamica EGO, calibrated with weights of evidence on eight variables describing mining and logging concessions, transport infrastructure, settlements and previous clearing. Forest covered 83.0% of the island in 2014 and 82.1% in 2024; under unchanged policy it falls to 73.7% by 2054, a net loss of 162 thousand ha, or 11.2% of the 2014 baseline, at gross rates of 47,000-51,000 ha per decade. Deforestation probability is highest within 500 m of previous clearing and declines with distance from settlements, cities and mining sites, while proximity to national parks carries a negative weight of evidence. The frontier is self-propagating and spatially predictable, and legally designated territory retains forest within it. Protecting the Hongana Manyawa consequently depends on excluding extractive licensing from the interior forest ahead of the frontier rather than behind it.

BackgroundVietnams reported dengue burden has risen roughly five-fold since 1990. Multi-decadal studies linking climate indices to dengue rarely separate genuine year-to-year coupling from a long-term trend the two share. MethodsWe assembled a provenance-preserving national annual dengue series (1990-2025; OpenDengue plus Ministry of Health figures) and correlated it with annual and March-May means of eight tropical sea surface temperature (SST) indices at lags of zero and one year under five trend-correction lenses: raw, linear detrending, first-differencing, socio-demographic-index residualisation and AR(1) prewhitening. FindingsCases rose 3{middle dot}5 percent annually. Seven predictors were significant at zero lag, led by the annual Indian Ocean Basin-Wide index (IOBW; Spearman +0{middle dot}534), but linear detrending removed all. First-differencing preserved six, led by spring IOBW (+0{middle dot}468), the annual Atlantic Multidecadal Oscillation (AMO; +0{middle dot}462) and annual IOBW (+0{middle dot}423); three survived AR(1) prewhitening - annual AMO and annual and spring IOBW. Spring AMO and a lag-1 Tropical North Atlantic signal (-0{middle dot}452) did not, and are hypothesis-generating. El Nino-Southern Oscillation indices failed throughout. InterpretationMost of the apparent association reflects a trend shared by warming oceans and expanding surveillance; we could not demonstrate that climate is the primary driver at this scale. Trend is not the whole story: IOBW and annual AMO persist under trend- and persistence-removing transformations. Because transmission responds to climate over weeks to months, annual averaging smooths the lags through which El Nino acts; these nulls reflect temporal scale, not climate insensitivity; usable predictors will require monthly, province-level models. FundingCenter for Environmental Intelligence, VinUniversity (project VUNI.CEI.FS_0001). Research in contextO_ST_ABSEvidence before this studyC_ST_ABSWe searched PubMed, Web of Science and Google Scholar for studies published up to May 2026 linking large-scale climate indices or sea surface temperature to dengue incidence, combining dengue, climate, sea surface temperature, ENSO, teleconnection and time-series terms with Vietnam, without language restriction. Many studies covering two or more decades reported strong correlations between basin-scale indices and national dengue counts. Most, however, relied on raw correlations or a single detrending choice, and rarely tested whether an apparent association reflected genuine year-to-year coupling or merely a shared long-term trend. Added value of this studyMost long-term studies remove the shared upward trend in only one way, or not at all. To our knowledge this is the first study to compare five trend-correction methods on a multi-decadal national dengue record and to read their agreement or disagreement as a diagnostic of which climate signals are real. A signal that appears only before the trend is removed is probably following it; one that persists is more likely real. Applied to a record spanning more than three decades, this comparison separates the two: several widely reported raw correlations weakened once the shared trend was accounted for. Implications of all the available evidenceClimate-informed analyses of multi-decadal data should report at least two trend-correction approaches alongside the raw correlation and treat their disagreement as evidence about where a signal sits, rather than operationalising raw long-span correlations. For Vietnam, the apparent national-scale association is dominated by a shared long-term trend but retains a smaller, robust inter-annual component led by the Indian Ocean and AMO signals; genuine coupling is more likely detectable at monthly resolution and provincial scale, where statistical power and physical mechanism are jointly available. Surveillance systems should retain explicit source provenance, so trend-corrected re-analysis remains possible as records grow.

Tropical cyclones are among the deadliest and costliest natural disasters in the United States, and the most intense storms are expected to become more frequent as the climate warms. Anticipating where deaths are most likely to occur is therefore central to preparedness, evacuation planning, and public health response. We modeled block-level mortality risk for twenty-four of the deadliest and costliest tropical cyclones to strike the U.S. Gulf and East Coasts, Puerto Rico, and the U.S. Virgin Islands between 1992 and 2024. For each storm, we combined NOAA hazard data (wind swaths, rainfall, and storm-surge inundation) with 2020 U.S. Census demographic and socioeconomic characteristics and the CDC/ATSDR Social Vulnerability Index for all Census blocks within 25 miles of the coast, and trained storm-specific boosted-tree models with population-standardized mortality as the outcome. Averaging block-level predictions within Saffir-Simpson categories yielded risk maps spanning tropical storms through Category 5 hurricanes. Predicted mortality risk rose with storm severity and concentrated in urban coastal communities of Puerto Rico, Louisiana, Florida, North Carolina, Virginia, Maryland, New Jersey, and New York, as well as in low-lying inlet, peninsula, and sound geographies. Large block population, non-Hispanic composition, male-dominated blocks, predominantly white blocks, and males aged 20 to 34 years ranked among the strongest predictors of mortality; patterns that likely reflect structural factors shaping exposure rather than individual susceptibility. The category-specific risk maps and an accompanying interactive dashboard provide a practical decision-support tool for emergency managers, planners, and coastal residents preparing for future storms.

Intensifying fire regimes threaten forests globally, but the risk of persistent post-fire forest loss and its potential mitigation remains poorly quantified. We analyzed millions of wildfires worldwide from 2001 to 2024 and tracked recovery in satellite-observed forest structure and ecosystem function. Post-fire persistent forest loss, indicated by modeled non-recovery to pre-fire conditions over decadal timescales, affected 57.1% of burned forest area globally since 2001, with hotspots in Pacific temperate and southern boreal forests. We then identified 'crucial fires' as events exceeding a stringent modeled-risk probability threshold for persistent structural or functional non-recovery, with fire severity strongly predicting this loss. This severity dependence revealed a management pathway, as locations with prior low-severity fire experienced lower severity in subsequent wildfires and had lower modeled probability of becoming crucial. Under a model-based counterfactual scenario, applying the estimated severity attenuation was associated with a 7.6% reduction; the top 1% of road-accessible areas accounted for 35% of this reduction. These results provide a global framework for identifying where wildfire threatens forest resistance and where targeted low-severity fire management like prescribed fire might be used to combat global forest loss.

Hydroclimate whiplash, defined as the rapid transition between unusually wet and dry conditions, is expected to intensify under climate change, yet its population health impacts remain largely unknown. Here we quantified the association between hydroclimate whiplash and mortality across the contiguous United States from 2003 to 2023 using monthly county-level mortality records, standardized precipitation evapotranspiration index data, and two-stage time-series models. We identified overall and direction-specific dry-to-wet and wet-to-dry whiplash events at seasonal and sub-annual timescales and across 5-, 10-, and 20-year recurrence intervals. More severe whiplash events were associated with higher all-cause mortality risk; 5-, 10-, and 20-year sub-annual overall whiplash events increased mortality risk over five months by 3.4%, 4.5%, and 5.7%, respectively. Elevated risks were observed across cause-specific mortality outcomes, with the strongest association for infectious diseases. We estimated that 103,471 deaths were attributable to overall whiplash during the study period. These findings identify hydroclimate whiplash as an emerging climate-related public health threat and suggest that adaptation strategies focused on single hazards may underestimate the health burden of rapid, sequential hydroclimatic extremes.

The tidal wetlands of Louisiana comprise about 25% of those found throughout the conterminous United States yet estimates of wetland loss rates in the region between 1932 and 2016 have exceeded 60 km2 yr-1. To mitigate further degradation and wetland loss in the region, a globally unprecedented $50B, 50-year plan for coastal Louisiana is driving restoration efforts, and demand exists from multiple stakeholders for regularly updated, regional-scale, accurate land cover information. We used machine learning (random forests; RF) and cloud computing to develop a new Landsat-based, marsh vegetation community geospatial dataset. The dataset depicts wetland vegetation community types defined in a previous study at annual (1985-2025) time steps at 30-m resolution. An RF algorithm was used to integrate training samples with feature variables derived from Landsat imagery, and the resulting geospatial data product achieved an overall correct classification rate of 78%. The approach for development of the land cover dataset presented here has potential for application in other coastal wetland habitats throughout the world.

Conserving migratory birds effectively requires full annual cycle strategies that identify where on-the-ground efforts can have the greatest impact. Here, we present a hemispheric spatial framework to identify priority areas for 112 migratory bird species across the Americas. Building on full annual cycle prioritizations, we defined finer-scale spatial planning units that reflect differences in migratory and congregational behaviors between shorebirds and landbirds. We compiled population data for each planning unit and focal species and applied conservation planning tools to design area-efficient portfolios of sites and landscapes that secure 10% of each species population within the Americas flyways. The resulting minimum area portfolios include 175 shorebird sites and 80 landbird landscapes optimized to meet the species-specific 10% representation targets across breeding, non-breeding, and passage seasons. We also identified a broader set of complementary solutions, ranked by an importance score, to provide decision-makers with flexible options for strategic resource allocation. This framework provides the scientific foundation for the Americas Flyways Initiative (AFI), which aims to catalyze investment in nature-based solutions and bird-friendly infrastructure to enhance the conservation of migratory birds and strengthen the resilience of the Americas flyways by 2050.

Invasive mosquito species pose substantial risks to human and animal health. Since 2004, Culex coronator, a mosquito vector species of public health concern, has shown rapid range expansion within the United States, spreading from a historically limited distribution in southern Texas to across the Gulf Coast region and into eastern and mid-Atlantic states. However, changes in environmental suitability associated with this expansion across historical, contemporary, and future climate conditions have not been evaluated. Here, we used species distribution models (SDMs) to compare predictions of abiotic suitability for Cx. coronator under historic (1960-1989) and recent (2000-2024) climate conditions calibrated on the historical range in the United States. We also created a contemporary SDM based on occurrence records prior to and following species range expansion (1960-2024), and further, to predict potential distributions under current and future climate conditions. Models calibrated on the historical range predicted only modest changes in suitability along the Gulf Coast region and failed to identify large areas of the humid subtropical eastern United States that are now occupied. In contrast, the contemporary model predicted widespread suitability across much of the southern and eastern United States. Future projections under the mid-range SSP3 scenario predicted increasing suitability at higher latitudes and elevations. Across all models, suitability was consistently low in arid and semi-arid regions, including along the historical western range limit, suggesting that moisture availability may constrain Cx. coronator distributions. Together, these results highlight the need to incorporate updated occurrence records when modeling invasive mosquito species to strengthen surveillance and control strategies.

Introduction Hearing loss (HL) affects over 1{middle dot}5 billion people globally and India shares a disproportionately high burden including Disabling Hearing Loss (DHL). HL affects an Individual socio-economically, but there are limited studies on the broader societal economic consequences of HL in India.Methods Using Cost-of-Illness (COI) approach, we studied the societal economic burden of HL in India. This study uses epidemiological and macroeconomic data and modelling to estimate the loss of Gross National Income (GNI) due to HL and DHL across three economic pathways. Uncertainty is evaluated using deterministic and Probabilistic Sensitivity Analyses (PSA).Results The model estimates that there are in India, 289 million and 85{middle dot}9 million people with HL and DHL respectively. Direct Loss of GNI and Indirect Loss of GNI (Caregiver burden) are estimated as INR 4,648{middle dot}4 billion (USD 55{middle dot}6 billion) and INR 3,268 billion (USD 39 billion) respectively. The Loss of GNI due to Low Education amongst those with HL is estimated as INR 1,041{middle dot}9 billion (USD 12{middle dot}45 billion).Discussion Economic burden of HL is presented across three pathways with Direct Loss of GNI due to DHL being the greatest. It also presents age stratified caregiver economic burden. The findings of the study help in estimating similar cost pathways, advocacy, and policy decisions towards reducing HL prevalence in India and LMICs. This study also highlights the need for India specific estimations related to the HL attributable low education, state-wise disaggregates, and prevalence studies. Funding This study has not received any funding.

Smallholder farming systems in sub-Saharan Africa are constrained by declining soil fertility, erosion, and rising fertilizer costs, creating an urgent need for scalable inputs that sustain yields while maintaining soil health. While there is some evidence that microbial inoculants may offer a promising complement to conventional fertility management, field-scale evidence in tropical cereal and tuber systems remains limited. Here, we evaluated a multi-species inoculant composed of 20-22 Bacillus and Streptomyces species on potato and maize across four sites in Rwanda over two growing seasons (2025A and 2025B). Treatments included the inoculant applied at two rates (150 and 250 g ha-1), both alone and in combination with standard fertilization (inorganic fertilizer plus manure), alongside untreated and fertilized controls. Co-application of the inoculant with standard fertilization increased yield and plant biomass beyond fertilization alone, with gains of 6-51% for maize and 3-58% for potato. However, while the inoculant applied alone outperformed untreated controls, it generally did not match standard fertilization. Responses were strongest and most consistent for large-grade potato tubers, and application rate interacted with crop type, whereby the lower dose maximized marketable tuber yield, while maize showed a positive dose-response for grain and biomass. Yield increases were not accompanied by reductions in crop nutrient density, which was instead governed by site-level differences. Altogether, these results indicate that multi-species microbial inoculants are an effective complement to existing fertility practices that may offer, pending further research, a potential pathway to partial fertilizer replacement while sustaining productivity and nutritional quality in smallholder tropical agriculture.

Background Childhood respiratory mortality in China has fallen by over 90% in three decades alongside sustained national warming, yet national long-run evidence on temperature and child respiratory mortality is lacking. Methods We linked Global Burden of Disease (GBD) 2021 mortality estimates for China - lower respiratory infections (LRI), ages 0-19, and asthma, ages 0-24, 1990-2021 - with C-LSAT 0.5 deg gridded temperature data (1990-2019), aggregated nationally and to five climate zones. Four annual indicators (mean temperature, diurnal temperature range, seasonal amplitude, interannual variability) entered regressions of log mortality rates with Newey-West standard errors. A bootstrapped (500 resamples) quadratic model probed the minimum mortality temperature (MMT), with PM2.5-adjusted analyses and future-exposure, permutation, and detrended falsification tests. Results LRI deaths fell by 96.3% (330,194 in 1990 to 12,098 in 2021; 95% uncertainty interval 9,669-14,891) and asthma deaths by 94.9% (3,287 to 167), while mean temperature rose 0.364 deg C per decade and diurnal temperature range narrowed 0.092 deg C per decade. Baseline coefficients were large (mean temperature -1.696, SE 0.174; diurnal temperature range +2.408, SE 0.336; seasonal amplitude -0.162, SE 0.082; interannual variability +2.924, SE 1.514, per 1 deg C in log rate), but the future-exposure test failed and detrending nullified every coefficient: the associations are trend-level, and short-cycle causal effects are not identifiable. Nor was the national MMT identifiable - observed temperature support spans only 6.66-8.13 deg C, and the nominal turning point of 35.84 deg C is an extrapolation artifact (quadratic term p = 0.963). Within the observed range, warming and declining mortality moved in the same direction. Conclusions The 96% decline in childhood respiratory mortality cannot be attributed to warming. China sits on the low-temperature side of the optimum, and the marginal direction of future warming requires stronger designs to establish. The falsification framework offers a discipline for climate-health inference in China.

Citizen science data are increasingly used to infer biodiversity change, but causal claims based on such data are credible only if sampling effort and its temporal shifts are explicitly modeled. Zhang et al. (1) used citizen science data to conclude that greater photovoltaic policy stringency, measured using the photovoltaic policy stringency index (PSI), reduced county-level bird diversity in China. We reproduced their fixed effects and instrumental variable estimates. However, the observed Shannon diversity derived from pooled citizen science records reflects both bird communities and sampling effort, which the authors' controls do not adequately capture. Accounting for observer count changed the reported statistically significant 2.10% decline in Shannon index to a nonsignificant 0.58% increase (P = 0.288) per one-standard-deviation increase in PSI, and rendered the instrumental variable estimate statistically indistinguishable from zero (P = 0.912). Yet observer count is only one of many sources of sampling bias. PSI was also associated with multiple dimensions of sampling effort, consistent with sampling effort acting as a potential mediator in the PSI-diversity chain. The sampling domain also shifted markedly from 2014 to 2023: recorded county-months increased almost 24-fold, median observer count rose from one to three, and zero-duration records declined from 57.2% to 0.17%. Without adequate adjustment, these shifts confound estimates of temporal change in observed bird diversity. Beyond its inadequate treatment of sampling effort, the original study also misinterpreted its statistical results. Although the reported R{superscript 2} values are high, they are dominated by county and year-month fixed effects, with PSI contributing a partial R{superscript 2} of only 0.048% on observed Shannon index. The PSI-photovoltaic-area correlation is also weak (r = 0.0414) and vanishes after accounting for fixed effects (P = 0.977). Furthermore, the released bird observation data contain many erroneous outliers, raising significant concerns about insufficiently rigorous data preprocessing and quality control. These results show that the released data cannot properly distinguish ecological change from sampling effort change. Robust inference from citizen science data requires checklist-level effort metadata, explicit correction for spatiotemporal sampling shifts, and close collaboration among researchers with complementary methodological and ecological expertise.

A direct empirical relationship between gross primary productivity (GPP) estimated by the eddy covariance method and satellite vegetation indices (VIs) has been widely observed across diverse ecosystems globally. Building on this observed covariation, VIs are frequently utilized as critical parameters - such as the fraction of absorbed photosynthetically active radiation (fPAR) - within light use efficiency (LUE) and greenness-based models for carbon cycle monitoring. However, actual canopy carbon assimilation is jointly governed by slowly evolving structural parameters and highly dynamic functional traits, such as physiological efficiency. The extent to which the macro-scale VI-GPP covariance is driven by structural scaffolding, and how this structural signal decouples from physiological function under environmental stress, remains to be systematically quantified. Here, we synthesized half-hourly eddy covariance measurements from 328 globally distributed sites and paired them with a rigorously angle-normalized Enhanced Vegetation Index (nadir view and fixed solar zenith angle at 30 degrees, EVI_SZA30). By applying a nonlinear light-response curve model across 54,720 high-frequency temporal windows, we mechanistically disentangled observed actual GPP (GPP_EC) into baseline photosynthetic capacity (P_c) and intrinsic quantum yield (alpha). Our results demonstrate that the macroscopic covariance between EVI_SZA30 and GPP_EC (R^2=0.554) is primarily driven by the index's robust ability to track structural capacity (P_c, R^2=0.538). In contrast, EVI_SZA30 exhibits limited sensitivity to high-frequency variations in functional traits like physiological efficiency (alpha, R^2=0.038). Particularly in water-limited biomes (e.g., open shrublands and woody savannas), intense environmental stress triggers rapid stomatal regulation while the physical canopy structure remains relatively stable. Consequently, the correlation between EVI and P_c becomes notably stronger than its correlation with actual GPP_EC, highlighting a pronounced structural-physiological decoupling. Because discrete overpasses by sun-synchronous polar-orbiting satellites face intrinsic temporal constraints in capturing sub-daily physiological down-regulation (e.g., midday photosynthetic depression), future monitoring paradigms could greatly benefit from the continuous, high-frequency observations provided by next-generation geostationary (GEO) satellites to bridge the gap between structural parameters and transient ecosystem function.

Many Sockeye salmon (Oncorhynchus nerka) populations have declined over recent decades, and climate change is likely to exacerbate these declines through direct and indirect ecological effects. The response to the associated environmental changes is likely to vary among life stages, populations, and regions. Quantitative estimates of climate change driven impacts that account for this variability could fill a critical gap and provide forward-looking insights into how sockeye are expected to respond to future climate-driven change across their lifecycle. To address this need we developed a hierarchical population dynamics model parameterized with juvenile, adult return and spawner abundance data from 13 sockeye salmon populations from Washington State to northern British Columbia. We used a formal causal inference framework that paired salmon abundance data with a suite of environmental covariates hypothesized to represent ecological conditions across the lifecycle. We used the model to estimate population-specific responses to each environmental driver, then combined parameter estimates with projections from down-scaled climate change models to estimate productivity responses to anticipated environmental change. We found that historical sockeye productivity was strongly associated with environmental covariates, which explained more interannual variability in return abundance than spawner abundance in most populations. However, the life stages and specific environmental covariates with the largest impacts differed among populations and regions, often displaying a latitudinal gradient. Increases in coastal ocean temperatures and mixed layer depth generally had negative effects though they varied among regions. Increased freshwater summer rearing and return migration temperatures had weaker but consistently negative effects. Under future climate conditions, projected changes in these environmental covariates are expected to result in substantial declines in productivity across most populations. Sockeye salmon display varying degrees of sensitivity to climate change across life stages, populations, and regions. Effective future management will require explicitly accounting for these life stage and population-specific responses.

Sequential hurricane-heatwave events threaten forest resilience via understudied legacy effects. Using Bayesian Structural Time Series, piecewise Structural Equation Modeling, and a 21-event global synthesis, we quantify how structural degradation non-linearly amplifies productivity loss during subsequent heatwaves. Our Hurricane Michael (2018) case study reveals significant negative GPP legacy effects during the 2019 heatwave. Intact, tall and diverse canopies buffer microclimates and moderate thermal sensitivity. Hurricane-induced structural simplification removes this protection, exposing temperature-sensitive shaded leaves to extreme stress. We identified a context-dependent hydraulic trade-off: structural complexity provides shading but exacerbates forest sensitivity to water deficits during peak heat, the vulnerability of which reverses during the recovery phase. Globally, these legacy effects are triggered by heatwave intensity and modulated by soil type, with loamy-soil forests most vulnerable. These findings highlight the critical role of forest structure in forest responses to compound disturbances. Neglecting structural legacies in Earth System Models likely underestimates risks to global carbon sinks.

High-resolution monitoring of forest structure and productivity is essential for effective natural resource management. However, monitoring approaches such as field-based forest inventories or extensive lidar campaigns are costly, time-intensive, and spatially limited. Therefore, inexpensive and accessible methods are needed. SatCHM (Satellite Canopy Height Model) was developed to be an accessible and open-source tool for researchers, allowing for site-specific and temporally flexible predictions of canopy height with limited computational resources. SatCHM requires four inputs: panchromatic satellite imagery, solar and sensor angle metadata of satellite imagery, digital elevation models (DEMs), and lidar-produced CHMs for an area of interest. After SatCHM pre-processes inputs, data is loaded into a collection of convolutional neural networks (CNNs) for image-to-image regression. This ensemble cooperates to yield high-resolution predictions (up to 0.5-meter) of three-dimensional tree structure with discernible tree crowns across a broader defined area of interest. After calculating the mean absolute error for each prediction output, the median of these mean absolute errors was 6.06 meters.